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Abstract 

Recent work has proposed the principle of 'asymptotic past triviality' to 
characterize the initial state in the pre-big bang scenario of string cosmology, 
that it is a generic perturbative solution of the low-energy effective action. 
Among the more generic sets of solutions which is simple enough to investigate 
thoroughly, yet complex enough to exhibit interesting behavior, is the gravity- 
dilaton system in spherical symmetry. Since, in the Einstein frame, this system 
reduces to a massless minimally coupled scalar, which has been target of a 
large body of previous investigation, we will draw on this and interpret it in 
the cosmological context. Since this scenario necessarily involves the transition 
from weak field initial data into the strong field regime, gravitational collapse, 
we have made numerical computations to answer some of the questions raised 
on the road to the proposal that 'the pre-big bang is as generic as gravitational 
collapse'. 
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1 Introduction 



Since superstring theory in its various manifestations JTJ, is currently the best can- 
didate for a theory incorporating General Relativity (GR) in a quantum context, it 
is natural to explore its implications for cosmology, in particular what it can say 
about a proposed early period of rapid expansion, inflation [[J. While string theory 
contains GR as a limiting case, the metric tensor is supplemented by other massless 
fields, universally including the scalar dilaton and an anti-symmetric tensor. There 
is additional field content dependent on the choice of a particular perturbative string 
theory, including numerous higher index form fields and moduli || . 

Of particular importance is the dilaton, which plays the role of a Brans-Dicke 
type field dynamically controlling the strength of the gravitational coupling. While 
the dilaton was recognized as a candidate inflaton early on, it was also realized that 
it created problems for standard potential dominated inflation Q. However it was 
also realized that a massless dilaton could drive inflationary expansion in the early 
universe (mass would be acquired later, presumably through a mechanism related 
to supersymmetry breaking) by virtue of an instability related solely to its kinetic 
energy. Using this as a base, Gabriele Veneziano and collaborators made an 'heretical' 
proposal for cosmological initial conditions 0, the pre-big bang (PBB). Since kinetic 
dilaton-driven inflation naturally pushed the universe in the direction of increasingly 
strong coupling and increasingly large curvature it was natural to propose that the 
universe began in a regime of weak coupling and small curvature, weakly coupled 
perturbative string theory in a nearly Minkowski vacuum. Furthermore, if a transition 
to a decelerated evolution and a stabilized coupling could be made, this era could be 
smoothly joined to a standard radiation dominated cosmology. While the mechanism 
of this transition ('graceful exit') remains obscure, it seems plausible that it could 
be achieved by high curvature and strong coupling corrections to the lowest order 
effective action ||, which are not presently unambiguously calculable. 

This radically different choice of initial conditions could not long escape the notice 
of cosmologists accustomed to thinking of initial conditions at near Planck scale 
conditions. In [0] it was pointed out that the initial spatial curvature scale of the 
universe would need to be extremely large in string units to give sufficient e-folds of 
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inflation, and || underlined this objection and quantified it in other ways, pointing 
out that it seemed to lead to other unnaturally large or small numbers and questioning 
whether the PBB in fact solved the problems inflation is supposed to solve. Part of 
these objections can be dealt with simply by pointing out that the posited PBB 
initial conditions constitute a classical state, so there is no reason why it should 
have a characteristic scale of the order of the string scale. Indeed the limit to the 
initial curvature scale is of the same order as its initial expansion rate ||, so there 
is simply no curvature problem to solve. Nonetheless, the question of how broad a 
range of initial conditions could be regarded as good for PBB inflation remained. In 



[Kj the authors explored the idea of a constant negative curvature Milne universe as 
a universal past attractor for PBB universes, but as pointed out in [inj a broadened 
class of Bianchi types includes solutions exhibiting gravitational plane waves as past 
attractors. 



This led to the proposal of |L2|, that a sufficiently generic state for PBB initial 
conditions is 'asymptotic past triviality' (APT), a generic perturbative solution to 
the lowest order effective action. More concretely, they propose a bath of weak 
gravitational and dilatonic radiation. They pointed out that this initial condition 
of radiation can still exhibit the runaway dilaton behavior characteristic of the PBB 
in the form of localized gravitational collapse. The interior of this collapsing region 
may serve as a seed (a baby universe). Just as in the simplest flat PBB, the runaway 
inflation produced by the kinetic energy of the dilaton is hoped, even expected, to be 
stopped by the singularity regulating effects expected of strings as extended objects, 
a graceful exit. With the abrupt curvature changes would come copious particle 
production and a transition to a decelerated evolution, a hot big-bang, followed by 
the standard cosmological scenario. 

A threat to this general notion of a black hole containing our universe is pioneer- 
ing work by Belinskii, Khalatnikov and Lifshitz (BKL) |TJ] suggesting that a generic 
approach to a vacuum singularity in Einstein gravity exhibits an infinite number of 
violent anisotropic oscillations. Encouraging is the observation by two of the same 



authors, among others [jlj], that the presence of a scalar field quenches these oscilla- 



tions. However, recent work |L5[ has suggested that truncating our field content to 
the graviton and dilaton is unrealistic, and that including the full form field content 



2 



of the full string theory will restore the BKL oscillatory nature of the singularity. 
A modified realization of this idea of realizing PBB initial conditions via collapse 



was proposed in [It]. They suggest a similar scenario in the context of singular- 
ity formation by colliding plane waves. By examining the analytic approach to the 
singularity cosmo logically in the string frame, they find a significant region in the 
parameter space of initial data which supports PBB inflation, in the sense of having 
all directions expanding in the string frame. But the effect of including the antisym- 
metric tensor field in the plane wave collision scenario was studied in ]17|], with the 
conclusion that this also tends to spoil the inflationary behavior and that PBB infla- 
tion may not be able to replace altogether more standard forms of potential driven 
inflation. Nevertheless, they suggest it may play a complementary role in providing 
initial conditions for subsequent standard inflation, so singularity onset from per- 
turbative initial conditions may still play an important part in the formation of the 
universe. 

Restricting the field content to the graviton and the dilaton and applying a con- 
formal rescaling from the string frame leads to a physically equivalent picture in the 
Einstein frame, gravity coupled minimally to a massless scalar. Further restriction 
to spherical symmetry gives a framework which is both theoretically and numerically 
tractable, but is still thought to retain the essence of physically realistic gravitational 
collapse. A long series of works by Christodoulou used this system to prove rig- 
orous results ranging from the existence of global solutions to Einstein's equations, 
to criteria guaranteeing the collapse of initial data to his most recent investigations 



of the genericity of naked singularities [lq|-|p3 . Further evidence of the richness of 
the system comes from the discovery of universal scaling and various sorts of critical 



phenomena 24 



We will concentrate on some questions raised in | l2"f . To probe the system the 
authors develop weak and strong field expansions for the evolution. However matching 
them only seems possible numerically, and we have developed an algorithm to do this 
matching. Our main goal is to explore the degree to which the black hole interior 
can resemble the dilaton driven regime of the PBB scenario. Other numerical studies 
of the PBB have focused on questions of its perturbative stability fl2~5| , but we will 
study the full nonlinear evolution. We will also look at a collapse criterion proposed 
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in [ 12 1 on the basis of a weak field expansion and attempt to draw some comparison 
with a more rigorous bound of Christodoulou. 



2 The Action, Conformal Frames and EOMs 

The tree-level low-energy effective action of a generic string theory, truncating to the 
graviton-dilaton content is 

S s = — ( d^xVGe-^lRiG) +G^d„$d u $] , (2.1) 
a J 

where G is the string frame metric, $ is the four dimensional dilaton and a' is related 
to the fundamental string length scale by a' ~ £ s 2 . While this is the physical frame 
in string theory and it is in this frame that the inflationary behavior of the dilaton is 
manifested most clearly, we will find it convenient to work in the conformally related 
but physically equivalent Einstein frame. Defining 

9^ = e^G^ (2.2) 

this action becomes 

Se = I [ d *x V# [R(g) - ~<T fy* d u $] , (2.3) 
a J 2 

the usual action for gravity coupled to a massless scalar. 

Reducing to spherical symmetry gives a numerically tractable but still rich sys- 
tem. We use an algorithm based on a form of the scalar-Einstein equations in radial 
symmetry given in a work of Christodoulou [|18| . Its use as the basis for a numerical 



code was pioneered in and has since been exploited by many authors [p7[] . It 
leads naturally to a completely characteristic stepping code. As such it has the same 
capabilities as double null codes [|28], , although it is most naturally expressed in a 



mixed coordinate system (u, r), where u is a null coordinate and r is the area r coor- 
dinate (spheres at fixed r have area Airr 2 ) . Though thought of as a coordinate in the 
definition of the metric, in the formulation of the algorithm r will become a quantity 
to be evolved along a second null direction. We deliberately do not specify whether u 
corresponds to an ingoing or outgoing null direction, since will have occasion to use 
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it both ways. Use of this algorithm along ingoing nulls was made in [31] to explore 
charged singularity structure. 

Through the introduction of an auxiliary variable we reduce the second order 
field equations to a first order system. Further, several of the field equations will 
be enforced by quadratures along nulls. Thus the algorithm requires no numerical 
differentiation during the evolution. This leads to an extremely stable 'pocket-sized' 
code, accounting for its popularity. Finally the auxiliary function will turn out to be 
proportional to the 'news' function at past null infinity, facilitating comparison with 

0- 

We express the metric in the (u, r) coordinate system as 

ds 2 = —g(u, r)g(u, r) du 2 — 2g(u, r) du dr + r 2 duj 2 . (2.4) 
In terms of the scalar normalization of [pl|, the Einstein-scalar equations are 



R P v = 8ttV m 0(m, r)V„0(u, r), (2.5) 
together with the scalar field equation 

VV^ = 0, (2.6) 

where is related to $ by 

$ = 4v^F0. (2.7) 
For notational reasons we will also write h = <fi. We then define the auxiliary variable 

(rh) r = h. (2.8) 



To connect this with the notation of T2], we recall that they define the 'news' function 
in the following way. The free field behavior of initial perturbative data is given by 

g = /QQ-/(«) (2J) 

r 

with u — t — t and v — t + r. The news is then defined by 

M(v) = /». (2.10) 
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Putting this together and correcting for normalization, we see that 

2y/nh — > J\f 

as u — > — oo at constant v (past null infinity or X" ). 
Finally we define the differential operator 

9 



(2-11) 



D u (f) = f, u - Ur 



(2-12) 

From the form of the metric ( |2.4| ), this is differentiation in the null direction comple- 
mentary to the nulls of constant u. That is, 



D u (f) 



du 



(2-13) 



where v is a complementary null coordinate to u. The component content of the 
equations of motion ( |2.5| ) and ( |2.6| ) can now be written very simply: 

(h-h) 2 



(log#),r 

{rg),r 
D u (h) 

= C 



47T- 



(g-g)(h-h) 



2r 



(9-9) 
2r 



log(^) 
9 



lirr(D u h)' 



(2.14) 
(2.15) 

(2.16) 
(2.17) 



As promised all derivatives are taken along null directions, so we now begin to think 
of all functions as functions of (u,v), though we will abuse notation somewhat and 
sometimes write them as functions of (u,r), which presents no problem since r(u,v) 
will be monotone in v. As usual, conservation imposes a redundancy on the equations 
of motion and one can be regarded as superfluous if we insure it is satisfied 'initially'. 
We choose to regard the most complicated equation ( |2.17| ) as the constraint and 
use the others to construct an algorithm. Some tedious work of differentiating and 
plugging in the other eoms shows, 

-2C 



(2.18) 



This has the solution £ = C/r 2 , implying that if the constraint is satisfied at any 
single point along the null, it will be satisfied everywhere along the null provided the 
other equations of motion are satisfied. 
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We choose an initial u = uq null. We specify initial data for r and h on a grid 
of points in v along the null. The choice of the v coordinate is completely arbitrary. 
From this data (|2.8| ) is used to construct h. We can then use Q2.14|) to construct 
g and Q2.15Q to construct g. Clearly initial conditions will be required somewhere 
along the null, but once this is done the constructions can carried out by simple 
quadrature. The choice of these initial conditions will be dictated by the need to 
solve the constraint. Once we have r, h, h, g, and g we can use ( [2.16| ) to step h to 
the next null at Uq + Am. We use the obvious evolution operator for r following from 

(EUP, 

D u (r) = =2 (2.19) 

to step the r values on the grid to the next null. Having new values of r and h on 
the next grid allows us to iterate the process. 

The location of the apparent horizon is given by the solution to the equation 

V M rW - - = 0. (2.20) 

Upon crossing this surface even the outgoing nulls become 'ingoing', in the sense of 
having a decreasing value of r along the null, signalling the formation of a trapped 
surface and inevitable collapse. Finally the Bondi mass function is given by 

m =£(!-£), (2.21) 
z 9 

so the apparent horizon can also be characterized by the condition 2m jr = 1. 



3 Stepping Modes 

This formulation of characteristic evolution, while technically simple, cannot cover the 
whole of the collapse spacetime in a single coordinate patch for several reasons. First, 
the metric components being evolved, g and g, are specific to an ingoing or outgoing 
null coordinate system. In an outgoing null coordinate system these components 
will diverge on approach to the apparent horizon, signalling a breakdown in the 
coordinate system. A global ingoing null system does not have this problem, but 
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is at best awkward to attach to r = 0, even before the formation of the r = 
singularity, due to mixed boundary conditions. Finally, on approach to X~ , r will 
diverge. The primary problem we wish to examine, the relationship of asymptotic 
news to asymptotic behavior near the singularity involves approaching all of these 
limits. 

To deal with this we split the spacetime into three regions and develop variants of 
the basic algorithm described in section 0. The first (w-mode) is one based on outgoing 
nulls and imposing boundary conditions at the regular r = origin. The second (v- 
mode) is based on ingoing nulls and takes boundary conditions on an outgoing null 
boundary from data extracted from the u-mode variant, allowing us to approach the 
r = singularity. Finally in the third (I~ -mode), we rescale the equations to put 
the X~ boundary at a finite position on the grid. Next we describe these modes in 
more detail. 



3.1 u-mode Stepping on Outgoing Nulls Touching r = 



Here we include r = as the initial point in the grid and identify u as a null coordinate 
which is constant along outgoing nulls, Fig. [I]. We will fix our initial data by specifying 
h(uo,v) and r(uo,v) along an initial outgoing null at finite u = u ('b' in Fig. [I]). 
Regularity at r = (V in Fig. |l|) requires that h(u,r = 0) = h(u,r = 0) (|2.8|) and 
g(u,r = 0) = g(u,r = 0) (|2.15| ). Further we will gauge fix u to be proper time at 
the origin r = 0, so g(u,r = 0) = g(u,r = 0) = 1. In light of (|2.18|) , we see that 
regularity at the origin implies £ = 0, so the constraint is satisfied everywhere as 
observed in fl8| . We note the explicit quadratures of ([2.8|), fl2.14| ) and ( |2.15|) subject 
to these boundary conditions, 



h(u, v) 



g{u,v) 



g(u,v) 



r(u, v)Jo 



r(u,v) 



exp 



1 



4tt 



h(u, r')dr' 
r M (h(u,r') - h(u,r')) 2 



dr' 



g(u, r')dr' . 



(3.1) 
(3.2) 
(3.3) 



r(w, v ) Jo 

The choice of v coordinate along the initial null is arbitrary and serves only to 
identify which grid point belongs to which ingoing null. In practice v is generally 
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r=0 



Figure 1: Evolution along outgoing nulls in a portion of a collapse spacetime. Shown 
are the singularity S, the apparent horizon AH and X~ . Initial conditions are fixed 
by regularity at the central axis v (r = 0) and initial data r and h are supplied along 
the outgoing nullb. The arrow indicates the direction of null stepping. Beginning at 
point a the stepped ingoing nulls are truncated near the apparent horizon. Initial data 
is then extracted along null e for use in the ingoing null version of the algorithm. 



taken from a X evolution (see section |3~3| ) so it corresponds to the usual v ~ t + r at 
past null infinity. The evolution of r and h to the next null is now carried out using 
( g7T§ and ( gTLj ). 



As we evolve, the r value of grid locations will decrease eventually becoming 
negative. The point then becomes unphysical and is dropped from the grid and a 
new r = point is interpolated. This is a consistent procedure since the quadratures 
( |3.3|) will maintain the correct regularity boundary conditions at the new origin. It is 
also reassuring to consider the limiting case of the evolution of extremely weak news 
data, the flat space wave equation. In this limit u = t — r, v = t + r and g — g — 1 so 
D u (h) = 0. Thus h{u,v) oc J\f{v) = f'iy). Substituting this into (|3.1|) and correcting 
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for normalizations gives 



$(w, v) = - M{v{u,r))dr 



-f V f'(v')dv' 
rJu 



f(v) ~ /(«) 



(3.4) 



r 



so we correctly reproduce the weak field evolution ( |2.9|) . In terms of our stepping 



algorithm this corresponds to our procedure of simply discarding points as they pass 
to negative values of r. 

We also monitor the value of 



along the grid. If \x approaches 1 (the point 'a' in Fig. [l|), then we are approaching 
the apparent horizon. This version of the algorithm cannot proceed past this point 
since r begins to decrease along an outgoing null as it crosses the apparent horizon, so 
the (u,r) coordinate system is not even one-to-one. The value of r at the grid point 
where this limit was exceeded is noted as an estimate of the final black hole size tbh, 
as is an estimate of the asymptotic Bondi mass, the value of m in the outermost grid 
point. 

We now proceed with the evolution by simply discarding the outer parts of the 
grid where a limit on \x is exceeded (notice \x = at the regular r = origin as can be 
seen from a power series expansion of the functions in r). This allows us to approach 
close to the appearance of the singularity at r = 0. The integration is halted when 
the outermost value of r is less than some r m j n and we extract initial data (u, r and h) 
along the ingoing null containing that point ('e' in Fig. P for use as initial conditions 
for the ingoing null phase of the evolution (see section p.2|) . 

3.2 v-mode Stepping on Ingoing Nulls Excluding r = 

In this mode we abandon our dependence on regularity at r = 0. Indeed, generally 
r = will be a singular point. The basic algorithm is identical with that described in 
section ^| but we now reinterpret the 'it' coordinate to be a conventional v coordinate 
constant along ingoing nulls. Our region of integration is now bounded by an ingoing 
null and an outgoing null, Fig. |[ We use the extracted data from the previous w-mode 
integration to provide initial data on the ingoing null ('e' in Fig. |]), u, r and h (which 



2m 




(3.5) 



r 
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is obtained from the extracted h by numerical differentiation). On the outgoing null 
boundary ('b' in Fig. |2|) we similarly extract values for r, h, g and g, but change 
g — > —g and g — > —g corresponding to the gauge change v — > —v. This gives us a v 
coordinate increasing along the outgoing null as implicitly defined by the analog of 



D v (r) = (3.1) 



Equations (|2.8| ), Q2.14Q and Q2.15Q are solved using the initial data on the outgoing 



null, instead of using the regularity conditions at r = 0. 

We now evolve just as before, but in the v direction, stepping the data on the 
ingoing initial null outwards, employing the same equations as in section Q with v 
substituted for u. These ensure the solution of all of Einstein's equations except for 
the constraint ( |2.17 ). But as we have argued through ( p. 18 ), we only need to show 



the constraint is satisfied at one point along each ingoing null. We choose this point 
to be its intersection with the outgoing null boundary. So we wish to show 

jg-g) | DM D v {g) 8nr(D v h) 2 _ Q 
2r g g g 

along the outgoing null, where 

D M = (3.3) 
By virtue of our definition of v coordinate ( |3.1| ) we can replace the D v operator in 

DM = ^(f) (3-4) 

where 

A-(/)=(f) M - (3.5) 

So D r (f) is the same as f >r in the (u, r) coordinate system. Now from ( |2.8|) , (|2.14|) 
and Q2.15| ) respectively (which were used to construct our original outgoing null data), 
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r=() 



Figure 2: Evolution along ingoing nulls in a collapse spacetime as in Fig. El. Initial 
data from the u-mode evolution is supplied along the ingoing null e and the outgoing 
null b. The arrow indicates the direction of null evolution. The ingoing nulls are 
truncated as they approach the r = singularity S. 

we have, 

^ ,tn h — h , 
D r (h) = (3.6) 

EM = J&z*>! (37) 

g r 

D r (g) = 9 —°-. (3.8) 
r 

We have written h instead of h to emphasize that it refers to the function h in the old 
coordinates (u,r). h and h are related to the invariant h = <p by r derivatives along 
different null directions. Substitution of these results into (|3~2l ) yields an identity. So 
we conclude we have consistent data for an integration of the eoms. 

According to a theorem of Christodoulou |2(J, the apparent horizon will extend 



outwards from the regular r = axis in an 'achronal' manner, spacelike except for 



12 



possible null parts. So provided we have penetrated deep enough in r with our ingoing 
initial null the evolving integration region should encounter the apparent horizon and 
pass inside. Notice that, unlike in the -u-mode scheme, g can pass through zero and 
change sign, g = signals the apparent horizon and there is no divergence in the 
metric functions in the f-mode gauge. Once an outgoing null crosses the horizon its 
r value begins to decrease as it falls in towards the singularity and we truncate the 
grid again at some fixed r min . 

Once the horizon appears we have grid points inside falling towards the singularity 
and those outside flying away, so the region around the horizon quickly empties. We 
replenish them using a 'point-splitting' scheme j2l| , simply interpolating new points 
to keep the grid populated. Also the u coordinate grid becomes extremely flat. To 
integrate to large values of v we would need a 'gauge-correction' and to increase 
efficiency a 'point-removal' scheme to eliminate unnecessary points. The extent of 
the region of spacetime we are interested in has not required these steps. 

3.3 T~ -mode Outgoing Null Evolution Excluding r = 

Finally, we introduce one more operating mode to bring past null infinity, X~ , into the 
grid, Fig. ^. We will rescale the auxiliary functions involved in the w-mode integration 
in such a way that derivatives and functions are well defined and finite at Z~ . To 
this end we introduce a compactified u coordinate and v related quantity, 

u = — — — (3.1) 

u 

v = u + 2r = 2r - \. (3.2) 

u 

v will be gauge fixed to equal to v at X~ {u = 0). This is, as usual, an arbitrary 
choice since it plays no essential role other than being a label, but ( |3.2| ) ensures that 
v approaches the usual sort of v ~ t + r coordinate at Z~ . We also restrict ourself to 
initial data such that J\f(v) — for v < Vq. This allows us to impose vacuum boundary 
conditions along the ingoing boundary null v = v (V in Fig. §) so we are not obliged 
to include the origin on the grid. To be clear we work in the (u, v) null coordinates, 
though again we occasionally abuse the notation and regard functions as depending 
on (u, v). v(u,v) will play a role similar to r in w-mode, it will be evolved. Initial 
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r=() 



Figure 3: Evolution along outgoing nulls from X in a collapse spacetime as in Fig. PI 
Initial data is given along a portion ofX~ , b, and vacuum initial data is set along 
the ingoing null v. The arrow indicates the direction of null stepping. Data at finite 
r is extracted along a final outgoing null e. 



data specified at I are v(0, v) and h(0,v). We now define the auxiliary functions, 
which are rescaled versions of (13. 31) modified to contain an inner cutoff radius, 



h(u, v) 
h(u, v) 



g(u,v) 



h(u, v 
uh = 



— 1 fv(u,v) _ 

h(u, v )dv 



g(u,v) = u(g-l) 



(1 + UV) Jv(u,v ) 
-1 



II 



u\g-l) 



rv(u,v) a(u ,V') +Uh(u 

exp Unu / . — —dv ) 

Jv(u,v ) (1 + V'U) 

I rv(u, v) 

1 ~l~ ~J\ J~J 



- 1 



g(u, v')dv 



(1 + UV) Jv(u,v ) 

In terms of these variables the evolution equations read 

Da = u 2 D u 



(3.3) 
(3.4) 

(3.5) 

(3.6) 



(3.7) 
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D t v = -tj (3.8) 

n = (39) 

1 +uv) 



X (u — 0) is now a regular null on the grid. Explicitly, at X ( |3.5|) can be replaced 
by the only non- vanishing term in its power series expansion in u, 

g(0,v) =-47rjJ5J )) /i(O J S0W. (3.10) 

The evolution is halted when the innermost point of the data reaches a specified 
finite r and data is extracted. It can then be connected through vacuum to the origin 
and evolution is continued in w-mode. 



4 Numerics, Tests and Verifications 

We use the lowest order algorithms still compatible with quadratic convergence. The 
r integrations along the nulls are carried out using the trapezoid rule, the null stepping 
is carried using a second order Runge-Kutta method and linear interpolation is used 
where necessary. Besides the advantage of simplicity, these methods work tolerably 
well even with data of low differentiability (e.g. |L2] (4.37)). 



It is useful to choose an analytic test case exhibiting features in common with 
scalar collapse. A simple case is spatially flat FRW scalar collapse, the simplest 
version of the PBB. While too simple to provide a good test of an (r, t) code, having 
no spatial gradients, when written in radial null coordinates we will see it has features 
in common with radial scalar collapse and has non-zero gradients in the null directions. 

We write the FRW metric 

ds 2 = -dt 2 + a(t) 2 (dR 2 + R 2 duj 2 ) (4.1) 
in terms of conformal time r), 

ds 2 = a(7]) 2 (-d7] 2 + dR 2 + R 2 du 2 ) (4.2) 

* - fwf- (4 - 3) 
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Introducing null coordinates denned by conformal time at the origin, u = r] — R and 
v = rj + R, we can write the PBB solution, 



<f>{v) 



a{rj) 




log(a(r?)), 



(4.4) 



(4.5) 



in terms of these null coordinates. Finally we replace the u coordinate with the u 
coordinate defined in terms of proper time at the origin 



where r]o = 3wo/2 is a chosen initial time (a(?7o) = 1)- Finally we change to the area 
r coordinate by R = r/a(u,v). 

The resulting solution has divergent news at X~ , requiring us to specify initial 
data at finite u, but the rest of the algorithm proceeds just as with a realistic collapse 
since this solution has a spacelike singularity at v = —u (rj = 0) hidden by a spacelike 
apparent horizon at v — —u/3, allowing us to test the horizon crossing and f-mode 
approach to the r = singularity. 

First we verify the order of convergence of the algorithm by checking it against the 
exact solution ( f4.5| ). We take initial r and h data from the exact solution at Uq = —2, 
—2 < v < 0.5. The data is then evolved to u — —1 and the results are compared with 
the exact solution. Fig. |] shows the logarithm (base 10) of the absolute value of the 
largest error, A, in <p on the grid versus the logarithm of the number of initial grid 
points n = 50, 100, 200, 400, 800. It is clearly almost exactly quadratically convergent. 

The full algorithm is then deployed to test the near r = approach to the singu- 
larity. Anticipating the asymptotic analysis of section [], we numerically extract the 
small r behavior of the dilaton near the singularity. As we will discuss this is expected 
to be of the form $(t>) = 2a(v) log(r), and since this is homogeneous and isotropic 
collapse, the value of a(v ) = — a/3 is expected. Using 400 points on the grid, the 
horizon first appears at r = 0.42 (marked by fx > 0.98) and the f-mode integration 
is pushed to r = 0.001. The results are shown in the second figure in Fig. f|. We 
show the values of a 2 (v ) extracted from the asymptotics of m and g. The small v 
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Figure 4: Left: the log of the maximum difference, (A), of the numerical integration of 
(j) from the exact PBB solution ( U-3j) plotted versus the log of the number of grid points 
n for n = 50, 100, 200, 400, 800. Right: a numerical extraction of the asymptotic value 
of a 2 (v) for the same solution from (f) (o), m (O) and g (x) over a range of v. 



deviation from the asymptotic values are expected, since the null geodesies have not 
yet penetrated close to the singularity, see Fig. @. In this particular run we find 
a(v) — > —1.728 compared with — Vo = —1.732, confirming we can reach near enough 
to the singularity to extract asymptotic values of a(v). As we will discuss further in 
section |^, this is the quantity which will characterize the cosmological nature of the 
collapse. 

Among other tests of the code, we repeat a calculation of the collapse of a initial 
pulse in and compare with the data presented in with good agreement. We 
note that although his integration proceeds several orders of magnitude deeper into 
the black hole than we do, he finds almost no change in the asymptotic a(v) over 
the last few orders of magnitude, giving us some confidence we can get a good idea 
of the asymptotic behavior with a relatively shallow penetration inside the apparent 



horizon. We also compare with other numerical results (e.g. [26]) and other exact 
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solutions. In addition we verify agreement between the various stepping modes on 
regions of overlap and numerically evaluate the constraint. All tests indicate the code 
produces stable integrations in agreement with others in the literature. 



5 Collapse Criteria 



In attacking the question of the structure of scalar collapse, the authors of |12| develop 
a perturbative expansion of the weak-field behavior in terms of powers of the strength 
of the news function, as defined by ( p. 10 ). They obtain a concise result for the 



quadratic approximation of the mass ratio [i = 2m/r for news data J\f(v) = f'(v), 

H2(u,v) = / ax j Ix) — . (5.1) 

v — u Ju [ v — u 

They further simplify by noting that if one defines the average of a function by 

1 f v 

(g)(u,v) = / g(x)dx, (5.2) 

V — U Ju 

and defines a variance by 

Var(g)(u,v) = ((g - (g)(u, v)) 2 )(u, v), (5.3) 

this becomes just 

fiv{u,v) = Var(jV)(u,v). (5.4) 

Since \i — 1 is the signal of the formation of an apparent horizon leading to 
inevitable collapse, this leads them to define the maximum of its quadratic approxi- 
mation 

BDV(J\f) = sup Var (AT) (u,v), (5.5) 

u,v 
u<.v 

and state that news data will collapse if it satisfies an inequality of the form 

BDV(JV) > C(jV), (5.6) 

where C(J\f) depends only on the shape of the news and not on its amplitude. So 
C(jV) measures one aspect of the departure from linearity of the evolution of J\f. The 
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usefulness of a criterion of this form is measured by the complexity of the functional 
C(N) and would be maximized if it could be replaced by a constant. But it could still 
be of some use if C(N) would be bounded for a family of 'reasonably generic' initial 
data and perhaps only depend weakly on Af. Aside from its simplicity the criterion is 
suggestive because it shares the exact scale invariance of the massless gravity-scalar 
system (r — > kr). The perturbative expansion leads them to estimate C(N) = 0(1). 
They also examine an exact collapse, [R3], and find C(Af) = 1/4 for this particular 
case. 

To examine this numerically, we select a family of news data depending on a 
number of parameters, map the boundary between collapse and dispersal in this 
parameter space and then evaluate BDV(M) along this boundary. Specifically we 
define a 'wavelet', 

n{v) = g C os( 27l{v - Vo)f )eM JV ~ V / ). (5-7) 
w w z 

Assigning J\f(v) = n(v) we see that this initial data has four parameters, amplitude a, 
frequency /, width w and center vq. vq is clearly irrelevant since translation of news 
data just shifts the value of u where the pulse approaches the origin. The irrelevance 
of w is less obvious, but comes from the scale invariance of the system. So shifting w 
can change the size of black hole formed, but cannot affect whether it will form. We 
have verified the numerics show this same insensitivity. So we will set these variables 
arbitrarily (vq = —2, w = 1). We then choose a value for / and locate a value 
of a where the data disperses and a larger value of a where it collapses. Searching 
by bisection we find the critical collapse amplitude ao(f) for this value of /. Doing 
this for a number of values of / allows us to collect a set of critical values of ao(/) 
for different values of /. Finally calculating BDV(M) for these critical collapse news 
functions gives us a cross-section of the shape of C(M) across this family of functions. 

Before doing this we check that the algorithm can accurately approach the collapse 
threshold by extracting the Choptuik exponent of the scaling of the black hole radius 
with the difference in amplitude from the critical amplitude 

r B H oc (a - a o y, (5.8) 

where the constant 7 is expected to depend only on the type of collapsing matter and 
not on the details of the waveform p4] . The results for / = news data are shown 
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Figure 5: Left: Choptuik scaling for f = 0, w = 1 news (M{v) = n'(v) ( \5. % )). 
Log of the black hole radius versus log of the difference between the amplitude a and 
the critical amplitude ao = 0.1896. A best fit line to the first five points is also 
shown with slope 0.373. Right: The C(N) functional appearing in the BDV criterion 
as a function of frequency f. On the upper trace M{v) oc n(v) and on the lower 
M(v) oc n'{v ). 



are shown in the first graph in Fig. || A best fit line to the first five points yields a 
slope of 7 = 0.373. This can be compared with other values in the literature, 0.376 
in the first reference in [24| or 0.374 in . 



We will use the wavelet data n(v) as initial data N(v). We also take a second 
family of news functions Af(v ) oc n'(v ) While looking substantially different for low 
frequencies where the derivative of n(v) is dominated by the gaussian component, 
at higher frequencies these functions will differ mainly by a phase coming from the 
oscillatory component. The dependence of C(Af) as a function of / for these two 
families is shown in the second graph in Fig. 

As one might expect, the phase difference at higher frequencies makes little differ- 
ence in collapse, as the numerics confirm by the convergence of the two curves. The 
most remarkable feature of the curves is their decline at high frequencies. If C(N) 



20 



1 ■ 

0.5 



1 1.5 





0.5 



1 



r u 

Figure 6: Collapse with M oc n'{v), f = 0, w = 1 with amplitude normalized to 
0.4- Left: the evolution of /i at several successive values of u approaching horizon 
formation. Right: The values of CHR(J\f,u) for these values ofu. 



has no lower bound, then in light of the scale invariance of the collapse this suggests 
that a very long wave train of very low amplitude news may in fact be gravitationally 
unstable to collapse. Unsystematic experiments have produced data which collapse 
at BDV(J\f) as low as 0.06, but numerical experiments alone cannot give a defini- 



tive answer. This complements analysis in [12fl suggesting that for extremely slowly 
decaying news data C(J\f) may go to zero. 

We also wish to make some comparison of this criteria with a rigorous sufficient 



criterion of Christodoulou p0 |. It states that if for some value of u 



2 Am 
Ar 



> 



— log 

r 2 



n 



2Ar 

where ri < r 2 ,r 2 < 3ri/2, Ar = r 2 — r\ and Am = m(u,r 2 ) — m(u,ri), then future 
collapse is inevitable. In physical terms it states that if the energy flux contained 
between two spheres is larger than a function involving the ratio of the two radii, 
then we can predict collapse. 



+ 6r i _ 1 

r 2 



(5.9) 
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Figure 7: Collapse with M oc n'(v), f — 0, w — 1 and amplitude normalized to 
0.2. Left: the evolution of fi at several successive values of u approaching horizon 
formation. Right: The values ofCHR(Af,u) for these values ofu. 



Defining Ci(Af,u,ri,r2) to be the left hand side of this inequality and 
C r (Af,u,ri,r 2 ) to be the right and 

CHR(M,u)= sup Ci{Af,u,r 1 ,r 2 )/C r (Af,u,r 1 ,r 2 ), (5.10) 

ri,r 2 
r 1 <r 2 <3r 1 /2 

we can then state the criterion as simply that CHR(Af,u) > 1 is sufficient to insure 
collapse. 

As this criterion needs to be applied at finite u (since m — * as u —>■ — oo, CHR 
can never predict collapse at large negative u) and the BDV criterion needs to be 
applied at X~ , a direct comparison is impossible. The best we can do to investigate 
its strength in predicting collapse is to ask how early it can predict collapse before 
the formation of the apparent horizon. 

To test this we form initial data Af oc n'(v ) for / = 0, w = 1 corresponding 
to a gaussian pulse in <p. Since the normalized amplitude for threshold collapse of 
this data set is found to be 0.1896, we first normalize the amplitude of M{y) to 0.40 
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(Fig. |6|) to investigate a heavy collapse and then to 0.20 (Fig. ^) to investigate a near 
critical collapse. In the first graph in each figure we show the profile of \i versus r 
for a number of values of u leading up to horizon formation. In the second graph 
we show the corresponding values of CHR(Af,u). In the first case we find that even 
for the heavier collapse the criterion predicts the appearance of the apparent horizon 
only rather near to its formation. In the second we do not get near enough to the 
horizon formation for the criterion to ever predict the collapse. This would seem to 
point to the strong nonlinearities involved in the horizon formation, particularly in 
the near critical region where a very small change in the news amplitude leads to a 
very large change in the largest value of \x eventually attained. 



6 Black Hole Interior Asymptotics and the PBB 

The proposal that 'the PBB is as generic as gravitational collapse' is the central pro- 
posal of |12|], and the measure of validity of this proposal lies in the interior of the 
black holes produced by massless scalar collapse. Since inside the black holes the r 
coordinate becomes timelike, the geometry there describes a spacelike hypersurface 
approaching a future singularity and the geometry is naturally interpreted as cos- 
mological. And since the dilaton, the scalar field, also generically shows runaway 
behavior, this suggests that the black hole interior when interpreted in the string 
frame may generically exhibit the pole-law inflation of the PBB scenario. With a 
graceful exit mechanism (possibly due to corrections to the lowest order perturbative 
action coming from large curvature and strong coupling) to regulate the approach to 
the singularity, the black hole interior may be regarded as a 'baby universe'. The first 
step is to check that the regime where the dynamics is dominated only by the lowest 
order action without corrections has some resemblance to the PBB scenario. 



Thanks to the work of it is known that the interior solutions for generic 
initial data (without additional field content) will be asymptotically velocity domi- 
nated, a 'quiescent cosmological singularity'. The final cosmological approach to the 
singularity will take the form of Kasner evolution, an FRW evolution modified by 
some degree of anisotropy. As the PBB form of inflation is known not to wash out 
anisotropy (as is clear in our case, since once we are into the Kasner regime the 
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degree of anisotropy freezes and remains constant), unless we take the opinion that 
the graceful exit will also cure this problem, we should look for features of generic 
gravitational collapse that will also select for isotropy. Working in spherical symme- 
try takes us part of the way to isotropy. since this forces two of the Kasner exponents 
to be equal, but this restriction is quite artificial. Nonetheless, if we can show col- 
lapse favors a third exponent equalling the other two, this would be some evidence a 
non-spherical collapse could also favor isotropy. 

In addition to relative isotropy, a 'good' PBB collapse should satisfy several re- 
quirements related to APT. The dilaton should be increasing with $ << and the 
spacetime curvature should be small at the entrance to the asymptotic regime to 
ensure the lowest order action is valid for a sufficiently large number of e-folds of dila- 
ton driven inflation. But we observe that any solution can be mapped to a solution 
of arbitrarily small curvature and weak coupling by applying a constant shift of the 
dilaton $ — > $ + C or performing a scale transformation $(r) — > $(&r) (changing 
the scalar curvature R — > R/k 2 ) without affecting the asymptotic Kasner behavior, 
as will be seen from the following analysis. 

In this scenario the need for a 'tuning' of these two classical moduli of the collapse 
is acknowledged. But due to the invariances of the system we should also expect 
collapses to occur on all size scales and couplings. Only sufficiently large and weakly 
coupled collapses will lead to universes such as we occupy. Others are regarded as 
stillborn [R5J. This also makes the question of the degree of fine-tuning subject to 
anthropic considerations, since the birth of a universe is no longer a unique event. 
With these invariances, the task left to a numeric code will be to analyze the relative 
degree of isotropy in the cosmology within the collapse as a function of the shape and 
amplitude of the initial news data. 

We first recall the small r asymptotics derived in ^] (and confirmed nu- 
merically in ||3U|| ) and rederive the relation between the asymptotics of the different 
functions. We start with the dominant terms in the growth of the scalar as we ap- 
proach the r = singularity, 

$(v, r) ~ A^kJi ~ 2a(v) logr + C^v). (6.1) 



24 



Using the denning relation for h (2.8) we easily have 



h(v,r)=h(v,r) + ^l. (6.2) 



Inserting this into ( 2.14j ) we find that to leading order, 



g(v,r) ~ C g (v)r a2 ^ (6.3) 

g(v,r) ~ (6.4) 
r 

Recalling the relation g/g = 1 — 2m/r, we have in the limit of small r, 

m(v,r) ~ C m (w)r" a2 ^. (6.5) 

So we see that nature of the evolution near the singularity is controlled entirely by 
az(v). We will next see how this a(v) is related to the Kasner exponents in the string 



and Einstein frame. We will follow the discussion of UI2| , filling in some details and 
omitting some more general discussions. 

Inserting the limiting behavior above into (|2.4j) we find 

ds 2 ~ r 01 *^- 1 dv 2 + r a2{v) dv dr + r 2 du 2 . (6.6) 

Here we have omitted all C(v) factors (in the asymptotic Kasner region they are 
locally nearly constant) and numerical constants since they can finally be absorbed 
into rescalings of the coordinates. By completing the square, this can be written as 

ds 2 ~ -^r a2{v)+1 dr 2 + (dv + r -dr) 2 r a ^ v) - x + r 2 doo 2 . (6.7) 

Changing to Einstein frame cosmic time via 

<x 2 (v)+3 

(-t B )~r^— , (6.8) 

this becomes (again omitting constants) 

ds E 2 ~ -dt E 2 + (-t E ) 2Xl d Xl 2 + (-t E ) 2X2 {dx 2 2 + dx 3 2 ) (6.9) 
$ ~ j(v)\og(-t E ), (6.10) 
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where X2 and £3 are the spacelike coordinates of spherical symmetry (in du 2 ) and x\ 
is the remaining spacelike direction and 

\(\ ft2 (^) ~ 1 WN 2 /X Mtj 

Ai(^) = 2 , v : Q , A 2 f = , 7 (u) = . 6.11 

cr(t>) + 3 cr(t>) + 3 cr(t>) + 3 

We complete the transition into the string frame by applying ( p.2|) and changing 
to string frame cosmic time by 

2 + 7(d) 

-*s = (-^) * (6.12) 



obtaining 



where 



ds s 2 ~ -dts 2 +(-t s ) 2ai da;i 2 + (-ts) 2aa (da:2 2 + dx3 2 ) (6.13) 
$ ~ tr(v)log(-t s ), (6.14) 



a» + 2a(w) - 1 . . 2(a(u) + i; 

«l(f) = 2(-.\ 1 o„Y„A I o > a2 ( W ) " 



a 2 (t») +2a(t») + 3' v/ a 2 {v) + 2a(u) + 3' 
a 2 (v) + 2a(v) + 3 

We can now state the requirement for a 'good' PBB collapse in terms of a(v). 
Consistent with inflationary behavior, we should at least require all directions are 
expanding (all string Kasner exponents negative) and an increasing dilaton. These 
requirements can be met by looking for regions where —1 — v^2 < a(v) < — 1 with 
a(v) = —V3 corresponding to exact isotropy. Since changing the sign of J\f(v) reverses 
the sign of a(v), we will be concerned with the magnitude of a(v). 

While it is clearly difficult to generalize about the infinite dimensional space of 
all possible initial data, we will present several data sets extracted from the many 
we have run looking for general trends. We choose a collapsing data set and extract 
values of a 2 {y) by fits to the small r data from (|6 . 1| ) , (|6.3|) and (|6.5|) . We will take 
agreement between these values to indicate that we are qualitatively in the asymptotic 
regime. 

In the first set, Fig. || we choose gaussian news, a = 0.35, / = 0, v = — 2 and w = 
0.2 (|5.7|) . This is somewhat over the critical amplitude of ao = 0.309. The first graph 
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Figure 8: Collapse of an amplitude 0.35 gaussian news pulse, left: h as a function of 
v atl^ , the end ofX~ evolution and near horizon formation (in order of increasing 
amplitude). Right: the extracted values of a 2 (v) from <f> (o), m (D) and g (x). 



shows the profile of h at X~ , at the end of the X~ evolution and at the point where 
fi > 0.99 and it shortly after forms a black hole of radius Tbh — 0.1349. We run the v- 
mode integration over the region near the r = singularity to extract the a 2 (v ) values 
from data going as deep as r min = 0.00008. As the second figure shows, this includes 
the region of significant deviation of a(v) from zero. At larger values of v, a(v) is 
running to its Schwarzschild value of zero, in accordance with the dictum 'black holes 
have no hair inside or out'. We do find a region of a 2 (v ) > 1, corresponding a minimal 
inflationary criterion of 'all-directions expanding', and further the maximum value of 
a 2 (v) ~ 3 reached corresponds to relative isotropy, although it is a narrow region 
very near to the initial v of the black hole formation. The innermost points represent 
ingoing nulls that fail to penetrate into the asymptotic region, so the apparent fall-off 
in the value of a(v) should not be taken seriously. We also note that the onset of 
black hole formation occurs as the end of the pulse is striking the origin. 

In the second set, Fig. || we raise the amplitude to a = 1.5, forming a black hole 
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Figure 9: Collapse of an amplitude 1.5 gaussian news pulse, left: h as a function of 
v atX" , the end ofX~ evolution and near horizon formation (in order of increasing 
amplitude). Right: the extracted values of a 2 (v) from <f> (o), m (D) and g (x). 
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of radius tbh = 3.495. We cut off at r min = 0.002. The qualitative results are 
quite similar to the first case, pointing again to a small relatively isotropic region 
corresponding to in a narrow region close to the onset of black hole formation. With 
the increased amplitude this now occurs when the leading surface of the pulse arrives 
at the origin. 

In the third set, Fig. [TO we extend a gaussian pulse of width 0.2 by adding a flat 



portion of v length 1.6 in the center of the pulse. We also tune the amplitude to 
0.33 so the apparent horizon forms as the flat part of the pulse is striking the origin, 
although the wave form is no longer flat by then. This time the closeness of the 
approach to a 2 = 3 is striking. 

Finally, in the fourth set we send in additional pulses after the first. We also vary 
the exact form of the initial news data to 



h(v) = asm ( ) 



w 



(6.16) 
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Figure 10: Collapse of a lengthened gaussian pulse, left: h as a function of v at X , 
the end ofX~ evolution and near horizon formation (in order of increasing amplitude) . 
Right: the extracted values of a 2 {v) from <f> (o), m (U) and g (x). 



with v = — 2, a = 0.4 and w — 1, and setting h — outside of four peaks of the news 
function. Again we find an initial approach to quasi-isotropy and notice that further 
incoming news data after the formation of the black hole does not have much effect 
on interior asymptotics. 

Comparison of these data with our criteria for a 'good' initial PBB leads to limited 
good news. Large regions of the black hole interior do not reach strong enough 
gradients in $ to meet the 'all directions expanding' criterion, much less isotropy. On 
the other hand there is a strong hint that there may be generically a small region 
very close the initial point of singularity formation where these requirements are met. 

Indeed we note that as we approach smaller and smaller r, closer and closer to the 
first appearance of the apparent horizon at the origin, the asymptotic a 2 (v) seems 
to approach 3 and seems to resist exceeding this value. This tempts us to conjecture 
that singularity onsets in spherical symmetry may approach exact isotropy. This is 
also suggested by the fact the regular r = origin is a point of exact isotropy before 
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Figure 11: Collapse of a wave train, left: h as a function of v at X~ , the end of 
T~ evolution and near horizon formation (in order of increasing amplitude) . Right: 
the extracted values of a 2 {v) from (j) (o), m (U) and g (x). 



the collapse, so this quality may extend into the singular region. While suggestive, we 
have been unable to verify this analytically. However, this would imply the occurrence 
of this near isotropic region may come from the exact spherical symmetry imposed, 
so it would not clearly extend to nonspherical collapse. 

This singularity onset region is also the most numerically difficult to reach and 
our results so far are hardly conclusive. Analyzing the exact solution (|12[ (4.37)) we 
find it corresponds to a 2 (v) = 1. But this solution is neither asymptotically flat nor 
differentiable at singularity onset, so it may not be representative of more physical 
collapse data. There also seems to be some trend for near critical collapses to reach 
higher levels of isotropy, but it is not clear whether this is because the isotropic regions 
are absent for heavier collapse or that we are simply unable to reach them. 

While we conclude the string frame near-isotropy region seems to be generically 
small (though in a way difficult to quantify exactly), it does seem to be present. 



Further, as suggested in | 3q| , an initially small isotropic region may grow preferentially 
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relative to its more anisotropic neighbors. Along a slice of constant string frame time 
(where tg = is the time of singularity, lacking a graceful exit), the volume of space 
per unit dv goes as 

Vol(t s ,v) = (-t s f a ^\ (6.17) 

This expansion rate is indeed maximized at a(v) = —a/3 with £aj(t>) = — v3. Since 
successful PBB inflation is expected to be accompanied by a large growth in $, in 
turn giving many e-folds in string time and in turn of the scale factor, this would seem 
to say that even an initially small region of isotropy could be exponentially enhanced 
over its neighbors having less negative exponents of growth rate. 

However, this process of comparing growth rates at different places at the same 
'time' requires singling out a particular foliation of spacelike hypersurfaces. For ex- 
ample, computing the string frame volume growth rates along slices of constant 
or constant r results in a preference for a(v ) with relatively strong anisotropy. Given 
this ambiguity we need to invoke some physical reason for a given choice of foliation. 
This choice should have something to do with the details of the final graceful exit 
mechanism but it is not clear to us how to make this more precise without choosing 
some family of observers to give the proper weight to different regions. Perhaps, fi- 
nally, these observers should correspond to the locations of future civilizations who 
could observe this anisotropy. For a Bayesian look at this possibility we refer again 
to [0, since it is not clear to us that the numerics can add quantitatively to this 



discussion. 

Nevertheless, we routinely find collapses featuring regions of a(v) < —1 and indeed 
approaching isotropy, confirming the general picture of 'ballooning patches of inflating 
space' in [O . 



7 Conclusions 

We have numerically investigated massless scalar collapse in spherical symmetry in 
the context of its possible role as a precursor for a PBB inflationary cosmology. 
We find evidence that such collapses do provide plausible precursors. Our most 
interesting conclusion based on this preliminary survey, that there seems to be a 
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region of substantial isotropy very near the initial appearance of the singularity on 
the central axis, will need to be investigated further. It lies in a region that is 
numerically difficult to reach and it would be much more satisfying to have some 
analytical understanding of this phenomenon, especially the degree to which it is 
particular to spherical symmetry. This conclusion also will need to be qualified to the 
degree to which the additional field content of string theory can spoil this behavior. 

We have also investigated a proposed collapse criterion |L2] . We find suggestions, 
again at the fringes of our numerical reach, that special forms of news data, specifically 
long trains of low amplitude waves, can collapse even though the criterion would 
suggest they should remain perturbative. In an extreme form this phenomenon would 
suggest that perhaps a very long train of arbitrarily low amplitude may be capable 
of collapse through a long slow buildup of non-linearity. While this seems somewhat 
physically unlikely it would be interesting to try to get some sort of analytical handle 
on this sort of behavior as well. 
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